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Summary 

A reverse transcriptase, polymerase chain reaction (RT- 
PCR) procedure was used to amplify a segment of the 
genome of turkey coronavirus (TCV) spanning portions 
of the matrix and nucleocapsid (MN) protein genes (ap- 
proximately 1.1 kb). The MN gene region of three epide- 
miologically distinct TCV strains (Minnesota, NC95, In- 
diana) was amplified, cloned into pUC19, and se- 
quenced. TCV MN gene sequences were compared with 
published sequences of other avian and mammalian co- 
ronaviruses. A high degree of similarity (>90%) was 
observed between the nucleotide, matrix protein, and 
nucleocapsid protein sequences of TCV strains and pub- 
lished sequences of infectious bronchitis virus (IBV). The 
matrix and nucleocapsid protein sequences of TCV had 
limited homology (<30%) with MN sequences of mam- 
malian coronaviruses. These results demonstrate a close 
genetic relationship between the avian coronaviruses, 
IBV and TCV. 
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Introduction 


The Coronaviridae are a large group of RNA-contain- 
ing viruses that infect a wide variety of avian and mam- 
malian species [1, 2]. The family is comprised of two gen- 
era, Coronavirus and Torovirus, that share similarities in 
morphology, genome organization and genome expres- 
sion [3]. The coronavirus genome consists of a positive 
sense, single-stranded RNA molecule that is 20-30 kb in 
size [3, 4]. Virions are enveloped, pleomorphic, 80-220 
nm in diameter, and have club-shaped surface projections 
approximately 20 nm in length [3, 4]. Three major struc- 
tural proteins are known, the surface (S) glycoprotein (90- 
180 kD), an intermembrane matrix (M) protein (20-35 
kD), and a nucleocapsid (N) protein (50-60 kD) [3, 4]. 

Coronaviruses initially were subdivided into four anti- 
genic groups based on differences determined primarily 
by immunofluorescent (FA), enzyme-linked immunosor- 
bent assay (ELISA), and immunoelectron microscopic 
(IEM) studies [1, 2, 5, 6]. Mammalian coronaviruses were 
shown to comprise antigenic groups | and 2, whereas avi- 
an coronaviruses, infectious bronchitis virus (IBV) and 
turkey coronavirus (TCV), comprised groups 3 and 4, 
respectively [1, 2, 5, 6]. IBV and TCV were determined to 
be antigenically distinct from each other and mammalian 
coronaviruses based on IEM, hemagglutination inhibition 
(HI) and virus-neutralization (VN) studies [7, 8]; how- 
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ever, subsequent antigenic and genomic analyses have 
questioned these taxonomic groupings. Studies done by 
Dea et al. [9] and later by Verbeek and Tijssen [10] indi- 
cated that TCV and bovine coronavirus (BCV), a mam- 
malian group 2 coronavirus, were closely related based on 
antigenic and nucleotide sequence analyses. In contrast to 
the findings of these studies, antigenic analyses by Guy et 
al. [11] indicated that the avian coronaviruses, IBV and 
TCV, are closely related and comprise a single antigenic 
group within the coronavirus genus [11]. 

TCV is the cause of an acute highly contagious enteric 
disease of turkeys referred to as bluecomb disease [12]. 
Bluecomb disease was first identified in turkeys in 1951 
and a coronavirus was identified as the cause of the dis- 
ease in 1973 [12]; however, antigenic and molecular char- 
acterization of TCV has been hampered by difficulties 
associated with in vitro culture of the virus. Minnesota 
strain, a reference TCV strain, and field isolates have 
been successfully propagated in embryonated turkey and 
embryonated chicken eggs by inoculation of the amniotic 
cavity [12]. Recently, Dea et al. [7] reported the adapta- 
tion and serial propagation of TCV in a human rectal ade- 
nocarcinoma (HRT) cell line. 

The present study was undertaken to examine relation- 
ships between TCV and other avian and mammalian co- 
ronaviruses based on nucleotide sequence analysis of 
TCV RNA. A well-characterized region of the coronavi- 
rus genome [13-15] comprised of nucleotides from both 
the matrix and nucleocapsid (MN) genes was used as a 
basis for comparison of three epidemiologically distinct 
TCV strains and selected coronaviruses. 


Materials and Methods 


Viruses 

TCV (Minnesota) was obtained from the American Type Culture 
Collection (Rockville, Md., USA). TCV (NC95) was isolated from 
turkeys in North Carolina [11]. TCV (Indiana) was isolated from tur- 
keys in Indiana and obtained from Tom Hooper, Purdue University, 
Dubois, Ind., USA. TCV strains were propagated by amniotic inocu- 
lation of embryonated turkey eggs [11]. IBV (Massachusetts) was 
obtained from SPAFAS, Inc. (Norwich, Conn., USA) and propagated 
in 9- to 11-day-old embryonated chicken eggs [16]. 


Virus Purification 

Template RNA for reverse transcriptase, polymerase chain reac- 
tion (RT-PCR) was obtained from partially purified preparations of 
TCV as follows. TCV-infected turkey embryo intestines were pre- 
pared as 10% (w/v) suspensions in TNE buffer (0.01 / Tris-hydro- 
chloride (pH 7.4), 0.1 44 NaCl, 1 mM EDTA), homogenized, and 
clarified by centrifugation at 8,000 g for 20 min. Polyethylene glycol 
(Sigma) was added to the supernatant fluid to a final concentration of 
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10%, then incubated overnight at 4 °C with gentle stirring. The pre- 
cipitate was collected by centrifugation at 8,000 g for 20 min. The 
pellet was resuspended in 1/10 of the original volume in TNE buffer 
pH 7.4, layered onto a 20-60% (w/w) sucrose step gradient, and cen- 
trifuged at 80,000 g for 2 h. The interphase was collected and viral 
particles were pelleted from the sucrose by centrifugation at 80,000 g 
for 1 h; the pellet was resuspended in RNase free water. IBV was 
similarly purified from virus-laden allantoic fluids. 

Template RNA also was obtained from TCV (NC95) which was 
purified by banding to density in sucrose gradients. TCV (NC95) was 
propagated and partially purified as described above. Material at the 
interface of the 20-60% (w/w) step gradient was resuspended in an 
equal volume TNE buffer pH 7.4 and layered onto a 20% to 60% 
(w/w) continuous sucrose gradient and centrifuged at 80,000 g over- 
night. The gradient was fractionated in 1-ml aliquots by upward dis- 
placement with a density gradient fractionator. Fractions were 
diluted in TNE buffer pH 7.4 and centrifuged at 80,000 g for 1.5 h; 
pellets were resuspended in TNE buffer pH 7.4. 


Nucleic Acid Purification 

Nucleic acid was harvested from purified virus by incubation in 
0.5% SDS (Gibco BRL) for 5 min at room temperature followed by a 
phenol-chloroform extraction was described [17]. Pellets were resus- 
pended in RNase free water. 


RT-PCR 

Viral cDNA was prepared using a commercially available RT 
reaction kit (Promega). The RT reaction consisted of RT buffer [100 
mM Tris-HCl (pH 8.8), 500 mM KCl, 1% Triton® X-100], 5 mv 
MgCl, 1 mM each deoxynucleoside triphosphate (dNTP), 20 units 
rRNasin, 12.5 units AMV RT, 0.5 ug random primers, 1 ug template 
RNA, and nuclease free water to 20 pl. The reaction was completed 
in a programmable thermal cycler at 42°C for 30 min, 99°C for 5 
min and 4°C for 5 min per the directions in the kit. 

PCR was conducted in the same tube on the resultant cDNA ina 
final volume of 100 yl. The PCR consisted of 2.5 units PWO DNA 
polymerase (Boehringer Mannheim), PCR buffer [100 mM Tris-HCl 
(pH 8.85), 250 mM KCl, 50 mM (NH4)SO4], 2.5 mM MgSO,, 100 ng 
of each custom primer EcoM and EcoN (Gibco BRL), and nuclease 
free water to volume. Samples were placed in a hot (94°C) thermal 
cycler and amplified as previously described [14] with the annealing 
temperature at 52°C for 2 min. PCR products were analyzed on a 
1% Ultra Pure Agarose (Gibco BRL) gel with 0.1 ug/ml ethidium 
bromide (Sigma). 

EcoM and EcoN primers were synthesized based on previously 
described primer sequences (MIBV, NIBV) for PCR amplification of 
IBV RNA [13, 14]. EcoM and EcoN primers were identical to MIBV 
and NIBV primers with the exception that restriction sites in these 
primers were replaced with EcoRI restriction sites. EcoM recognizes 
the 3’ end of the matrix protein gene and EcoN recognizes the 5’ end 
of the nucleocapsid protein gene. Primer sequences are as follows: 

EcoM 5’-TGAATTCTCAGTGGCTTGCTAAGTGTGAACC-3’ 

EcoN 5’-TGAATTCACCGCTACCTTCAAACTTGGGCGG-3’ 


Cloning 

RT-PCR products and cloning vector pUC19 were digested with 
restriction endonuclease EcoRI (Gibco BRL) according to manufac- 
turer’s recommendations. RT-PCR products were inserted into 
pUC19 and transformed into a competent Escherichia coli-strain 
DH5a, (Gibco BRL) as described [17]. 
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Sequencing/Phylogenetic Tree Analysis strains), BCV (Mebus, F15 strains), mouse hepatitis virus (MHV) 


DNA was sequenced at the University of North Carolina(Chapel (A59, JHM strains), transmissible gastroenteritis virus (TGEV) (Pur- 
Hill) Automated DNA Sequencing Facility, on a Model 373A DNA due, FS772 strains), porcine respiratory coronavirus (PRCV) (86- 
Sequencer (Applied Biosystems) using the Taq DyeDeoxy® Termi- 137004, RM4 strains), canine coronavirus (CCV) (Insave-1 strain), 
nator Cycle Sequencing Kit (Applied Biosystems). All sequences human coronavirus (229E, OC43 strains), feline infectious peritoni- 
were confirmed by sequencing both strands. GenBank accession — tis virus (79-1146), and feline enteric coronavirus (79-1683) [15, 
numbers are: TCV (Minnesota) AF072911; TCV (NC95) AF072912; 19-34]. The sequence reported by Verbeek and Tijsen for TCV, and 
TCV (Indiana) AF072913. herein referred to as TCV (Verbeek), was also included in the study 

Comparative analyses of nucleotide and protein sequences were [10]. 
performed using the MegAlign application of the Lasergene software 
package (DNASTAR, Madison, Wisc., USA). Phylogenetic tree con- 
struction was based on the neighbor-joining method of Thompson et 


al. [18] using an unrooted tree analysis in the program CLUSTAL X Results 
with 1,000 bootstrap trials. Bootstrap values =950 are considered 
significant at the 95% confidence level (p <0.05). Protein compari- RT-PCR 


son consisted of 70 amino acids at the carboxy-terminus of the a one OR ; 
matrix gene and 55 amino acids at the amino-terminus of the nucleo- Preliminary studies indicated that RNA obtained from 


capsid gene. Coronavirus sequences were obtained from the Gen- turkey embryo-propagated TCV (NC95) could be ampli- 
Bank database. These included IBV (Beaudette, Gray, KB8523 fied in an RT-PCR using IBV-specific synthetic primers 
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IBV gt [20] and TCV strains NC95, Minnesota, and Indiana. 

NC95 .. Nucleotide sequence differences are shown for TCV 

sa ae strains. The positions where nucleotide bases are miss- 


ing are indicated as (—) and similar nucleotides as (.). 


MIBV and NIBV (data not shown).These primers were 
previously described by Andreason et al. [13] for RT-PCR 


described by Andreason et al. [13]. No PCR product was 
observed when RNA was harvested from uninfected em- 


amplification of the MN gene region of IBV RNA. Subse- 
quent RT-PCR of TCV strains utilized synthetic DNA 
primers EcoM and EcoN; these primers were identical to 
MIBV and NIBV with the exception that EcoRI restric- 
tion sites were substituted for restriction sites in these 
primers. 

RNA obtained from partially purified TCV (NC95), 
TCV (Indiana) and TCV (Minnesota) was used as tem- 
plate in RT-PCR using synthetic primers EcoM and 
EcoN. The PCR products were approximately 1.1 kb in 
size and slightly larger than the product obtained using 
IBV RNA as template (fig. 1). IBV RNA was used as a 
positive control in the RT-PCR procedure and produced 
a product of approximately 1.0 kb, consistent with that 


Sequence Analysis of Turkey Coronavirus 


bryos, and amplified by RT-PCR, or when RT-PCR was 
run without RT (data not shown). 


RT-PCR Analysis of Sucrose-Gradient-Purified TCV 

As an additional means of validating the RT-PCR, 
template RNA was obtained from sucrose-gradient-puri- 
fied TCV (NC95). RT-PCR was performed on RNA 
extracted from individual gradient fractions after isopyc- 
nic centrifugation of TCV (NC95). Positive PCR amplifi- 
cation was evident in fractions with densities of 1.14-1.20 
g/cm}, consistent with the density of coronaviruses [3]. 
Positive PCR amplification was seen in only one other 
fraction, at a density of 1.113 g/cm?. 
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for TCV strains. The positions where amino 
acids are similar are indicated as (.). 


Fig. 4. Phylogenetic relationship of avian 
and mammalian coronaviruses based on 
comparison of matrix protein (A) and nu- 
cleocapsid protein sequences (B) in the 1.1- 
kb MN gene region. Amino acid sequences 
were aligned using the CLUSTAL X method 
[18] and unrooted phylogenetic trees were 
constructed using the neighbor-joining 
method with 1,000 bootstrap trials. The 
bootstrap values are indicated at the 17 in- 
ternal branch points. 


Sequencing/Phylogenetic Trees 

EcoRI digestion of the TCV RT-PCR product, in prep- 
aration for cloning, revealed an interior EcoRI site in 
both TCV (NC95) and TCV (Indiana), but not in TCV 
(Minnesota). Therefore, TCV (NC95) and TCV (Indiana) 
were cloned and sequenced in two pieces, one approxi- 
mately 400 bases and the other approximately 700 bases. 
TCV (Minnesota) was cloned and sequenced as one 1.1 kb 
piece. Figure 2 compares the MN gene nucleotide se- 
quence (1.1 kb) of three strains of TCV with the published 
sequence of IBV (Beaudette) [20]. Figure 3 compares the 
matrix and nucleocapsid amino acid sequence of three 
strains of TCV with the published sequence of IBV (Beau- 
dette) [20]. 

MN nucleotide sequences of TCV strains are com- 
pared in table 1 with published sequences of IBV strains 
(Beaudette, KB8523) and two representative mammalian 
coronaviruses, BCV (Mebus) and TGEV (Purdue); results 
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are shown as the percentage of sequence homology be- 
tween viruses. Nuclotide sequence homology among the 
three TCV ranged from 87.2 to 89.6%. Nucleotide se- 
quence homology between TCV strains and IBV strains 
ranged from 85.8 to 92.6%. Homology between TCV 
strains and mammalian coronaviruses was less than 
30%. 

Comparisons between matrix and nucleocapsid pro- 
tein sequences are shown in table 2. Matrix and nucleo- 
capsid protein homologies between TCV and IBV ranged 
from 95.6 to 100% and 90.9 to 98.2%, respectively. In 
contrast, homology between TCV matrix and nucleocap- 
sid protein sequences and those of mammalian coronavi- 
ruses was less than 30%. Similarly, homology between 
matrix and nucleocapsid proteins of TCV strains exam- 
ined in this study (Minnesota, NC95, Indiana) and TCV 
(Verbeek) was less than 30%. These findings demonstrate 
that TCV (Minnesota, NC95, Indiana) and IBV are highly 
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Table 1. Percent nucleotide sequence homology in MN gene region between TCV (Minnesota, NC95, Indiana), IBV 
(Beaudette, KB8523) [20, 21], and representative mammalian coronaviruses, TGEV (Purdue) [27, 28] and BCV 


(Mebus) [15] 
IBV(Beaudette) IBV(KB8523) TCV(NC95) TCV(Minnesota) TCV(Indiana) _BCV(Mebus) 
IBV(KB8523) 90.3 
TCV(NC95) 87.0 90.4 
TCV(Minnesota) 90.0 92.6 88.5 
TCV(ndiana) 85.8 88.4 89.6 87.2 
BCV(Mebus) 21.6 21.5 21.0 21.8 19.6 
TGEV(Purdue) 21.1 22.7 20.9 20.5 21.2 27.5 


Nucleotide sequence information was obtained from the GenBank database. Comparative analyses of nucleotide 
sequences were performed using the MegAlign application of the Lasergene software package (DNASTAR, Madison, 
Wisc.). 


Table 2. Amino acid sequence comparison of TCV matrix and nucleocapsid proteins with corresponding regions of 
other avian and mammalian coronaviruses 


Nucleocapsid percent similarity 


IBV IBV TCV TCV TCV TCV BCV MHV CCV TGEV 
(Beaudette) (KB8523) (NC95) (Minnesota) (Indiana) (Verbeek) (Mebus) (A-59) (Insave-1) (Purdue) 
IBV 93.4 90.9 96.4 90.9 20.3 20.0 21.3 19.7 20.4 
(Beaudette) 
IBV 95.6 92.7 98.2 92.7 20.3 18.3 21.3 19.7 20.2 
(KB8523) 
TCV 100.0 98.6 94.5 92.7 20.0 20.0 21.8 18.2 20.0 
(NC95) 
TCV 98.6 97.1 98.6 94.5 20.0 20.0 20.0 18.2 20.0 
(Minnesota) 
TCV 98.6 97.1 98.6 97.1 20.0 20.0 18.2 18.2 21.8 
(Indiana) 
TCV 28.4 27.6 19.6 17.6 19.6 99.6 67.4 17.1 19.9 
(Verbeek) 
BCV 28.4 28.0 19.6 17.6 19.6 99.6 67.2 17.1 19.9 
(Mebus) 
MHV 28.9 28.4 17.6 17.6 17.6 83.8 84.2 21,5 23.0 
(A-59) 
CCV 19.6 19.1 23.5 23.5 23.5 34.3 34.8 32.0 89.0 
(Insave-1) 
TGEV 20.4 19.6 23.5 23.5 23.5 36.5 37.0 35.1 86.3 
(Purdue) 


Matrix percent similarity. 


related (>85% at the nucleotide level; >90% at the level of prepared based on matrix and nucleocapsid protein se- 


the matrix and nucleocapsid proteins). 

Phylogenetic trees were prepared to examine relation- 
ships between TCV and other coronaviruses based on 
comparison of amino acid sequences within the matrix 
and nucleocapsid proteins (fig. 4). Phylogenetic trees were 
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quences for TCV strains and published sequence data for 
selected avian and mammalian coronaviruses. These phy- 
logenetic trees divided the coronaviruses into three geno- 
typic groups with the avian coronaviruses, IBV and TCV, 
representing one distinct group. 
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Discussion 


Previous studies in our laboratory indicated a close 
antigenic relationship between IBV and TCV [11]. TCV 
(NC95) and TCV (Minnesota) were shown to be closely 
related to IBV based on FA procedures using polyclonal 
and monolconal antibodies [11]. These studies suggested 
that these avian coronaviruses comprise a single antigenic 
group within the Coronavirus genus. The present study 
provides additional support for those findings. Sequenc- 
ing data demonstrated a close genetic relationship be- 
tween IBV and TCV, and phylogenetic analyses divided 
the avian and mammalian coronaviruses into three dis- 
tinct genotypic groups with the avian coronaviruses, IBV 
and TCV, comprising a distinct genotype. 

Synthetic primers that previously were described for 
RT-PCR amplification of IBV RNA within the MN gene 
region were shown in this study to amplify TCV RNA 
[13]. The RT-PCR product was very similar in size (1,100 
bp) to the expected product; RT-PCR amplification of 
IBV RNA resulted in a 1,000-bp DNA product [13]. Sub- 
sequent nucleotide sequence analyses of the resultant 
DNA products demonstrated that the sequence of the MN 
gene region of TCV was closely related to IBV, and only 
distantly related to mammalian coronaviruses. A high 
degree of similarity (>90%) was observed between the 
matrix protein and nucleocapsid protein sequences of 
TCV strains and sequences of IBV strains. In contrast, a 
significantly lower degree of homology (<30) was evident 
between matrix and nucleocapsid protein sequences of 
TCV and mammalian coronaviruses. 

Several control procedures were used to validate the 
RT-PCR procedure for amplification of TCV RNA. No 
RT-PCR product was obtained when RT-PCR was per- 
formed with material obtained from uninfected embryo 
intestines, and no product was obtained when the RT- 
PCR was performed without RT. The source of template 
RNA was examined by isopycnic centrifugation of TCV 
(NC95) in a sucrose gradient, followed by gradient frac- 
tionation, nucleic acid extraction, and RT-PCR amplifi- 
cation of individual gradient fractions. Gradient frac- 
tions which yielded positive RT-PCR products had den- 
sities of 1.113 g/cm} and 1.14 to 1.20 g/cm3. Coronavi- 
ruses have a buoyant density in sucrose of approximately 
1.15-1.19 g/cm? [3]. Thus, the material used for RNA 
extraction and RT-PCR banded at densities consistent 
with that of coronaviruses. The nature of the material 
that banded at a density of 1.113 g/cm? and produced a 
positive RT-PCR product has not been determined; 
however, this material likely resulted from virion degra- 
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dation and release of viral RNA during sedimentation in 
sucrose. 

Studies conducted by Dea et al. [9] and Verbeek and 
Tissen [10] indicated a close antigenic and genomic rela- 
tionship between TCV and BCV; however, our studies 
contradict those findings. In the studies by Dea et al. [9], 
TCV was shown to replicate in HRT cells and to be closely 
related to BCV based on immunoblotting studies, HI, 
VN, and hybridization of BCV cDNA to TCV RNA. Ver- 
beek and Tijssen [10] also used HRT-propagated TCV 
strains and their studies demonstrated nucleotide se- 
quence homology of 99% between TCV and BCV in the 
MN gene region. In contrast to the findings of Dea et al. 
[9] and Verbeek and Tijssen [10], our previous studies 
failed to detect antigenic similarity between BCV and 
TCV using cross-immunofluorescence, and our attempts 
to propagate TCV in HRT cells were unsuccessful [11]. In 
the present studies, sequence analyses supported our pre- 
vious antigenic analyses and indicated that TCV strains 
NC95, Minnesota, and Indiana were closely related to 
IBV. 

The results of the present studies demonstrate a strong 
genetic relationship between the avian coronaviruses, 
IBV and TCV, and support the contention that IBV and 
TCV comprise a single antigenic group/genotype within 
the Coronavirus genus. Further studies examining the 
genetic relationship between the surface glycoprotein (S) 
genes of TCV and other coronaviruses are warranted. 
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